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ABSTRACT 

We present Very Large Array (VLA) observations of the water maser emission 
towards IRAS 16552-3050. The maser emission shows a velocity spread of ^ 
170 km s~^, and a bipolar distribution with a separation between the red and 
blueshifted groups of ~ 0.08". These observations and the likely post-AGB nature 
of the source indicate that IRAS 16552-3050 can be considered as a member of 
the "water fountain" class of sources (evolved stars showing H2O maser emission 
with a velocity spread > 100 km s~^, probably tracing collimated jets). The 
water maser emission in IRAS 16552-3050 does not seem to be associated with 
with any known optical counterpart. Moreover, this source does not have a near- 
IR 2MASS counterpart, as it happens in about half of the water fountains known. 
This suggests that these sources tend to be heavily obscured objects, probably 
with massive precursors (^4 — 5 M©). We suggest that the water maser emission 
in IRAS 16552-3050 could be tracing a rapidly precessing bipolar jet. 

Subject headings: masers - stars: AGB and post-AGB - stars: individual (IRAS 
16552—3050) - stars: mass loss - stars: winds, outflows 



1. Introduction 



Water fountains are evolved Asymptotic Giant Branch (AGB) and young post-AGB 
stars that sho w water maser emissi on with velocity separations in their spectral features 
>100 km s~^ ( iLikkel &: Morrislll988l ). Some of them have been observed with radio interfer- 
ometric techniques and all show bipolar distributions of their H2O mas er emission, indicating 



the presence of jets with extremely short dynamical ages of 5-100 yr (ILikkel fc Morris! Il988 



Boboltz fc Marvell l2005l . 120071 : llmai et al.ll200l 12004 : llmail 120071 ). OH maser emission has 
been detected in most sources catalogued as water fountains and, in general, this emission 
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exhibits lower velo cities and is spatially less extended than that of H2O (see, for example 



Deacon et al.l 120071 ) . Water fountains are evolutionarily located between two phases of stel- 



lar evolution, AGB stars and planetary nebulae (PNe), when mass- loss properties, among 
others, change dramatically. While mass ejection during the AGB phase is spherical, most 
PNe exhibit elliptical, bipolar, or multipolar morphologies. These morphologies have been 
attrib uted to the action of c oUimated outflows on the previously expelled spherical AGB 
shell ( ISahai fc Traugeiill998l ). Therefore, water fountains very probably trace the onset of 
non-spherical and highly collimated mass ejection. Understanding the transformation of an 
AGB star into a PN relies on a precise knowledge of the properties and evolution of water 
fountains. Unfortunat ely, very few water fountains are known at present (11 sources reported 
as of March 2008, see llmail 120071 ). so that adding new members to this important class of 
objects represents a valuable progress. 



Recently, in a single-dish search for H2O masers in evolved stars, ISuarez et al.l (120071 ) 
found water maser emission toward IRAS 16552—3050 (hereafter IRAS16552), with a maxi- 
mum separation among its velocity comp onents of ^170 km s~^ . IRAS 16552 was first pro- 



posed to be a candidate post- AGB star by IPreite-MartinezI (1l988l ). based on its IR AS colors. 



A pos sible optical counterpart of the IRAS far-infrared source was classified by IHu et al. 



( I1993I ) and ISuarez et al.l (120061 ) as a star of spectral type K9III and KOI, respectively. 



In this paper, we present Very Large Array (VLA) observations, carried out to confirm 
the association of this high-velocity water maser emission with the IRAS source, and to 
determine the spatial distribution of the masers. Our results allow us to include IRAS 16552 
as a bona fide member of the water fountain class. In Sec. [2] we describe the performed 
observations. Sec. [3] presents the obtained results, and Sec. |l]a discussion on the properties 
of IRAS16552 and a comparison of this source with other water fountains. Finally, in Sec. 
E]we summarize the main conclusions of this work. 



2. Observations and data reduction 

We have observed the 616 — 523 maser transition of the water molecule (rest frequency 
= 22235.08 MHz) using the VLA of the National Radio Astronomy Observatory (NRAO)0 
in two separate epochs. In the first one (2006 March 22 and 26), we used the A configuration 
of the VLA, and observed the right circular polarization (1 IF mode) over a bandwidth of 
3.11 MHz sampled with 255 channels (thus providing a spectral resolution of 1.22 kHz, or 
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0.1 6 km s To be a ble to cover the whole range of the water maser emission detected 



by ISuarez et al.l (120071 ) toward IRAS16552, we successively centered the bandwidth at six 
different velocities with respect to the Local Standard of Rest (LSR): —70, —30, and +10 
km s-i on 2006 March 22, and +50, +85, and +120 km s'^ on 2006 March 26. Therefore, 
we covered the water maser emission of IRAS16552 on a velocity range from —91 to +141 
km s~^ (although not simultaneously). The phase center for these observations was set 
at R.A.(J2000) = 16''58'^27.84^^ecXJ2000}.^-30°55'06.7'', which is within 1" from the 



position of the source given in ISuarez et al.l (120061 ) . 



In the second epoch, carried out on 2007 December 24, we used the B configuration of the 
VLA and the 4 IF mode, with a bandwidth of 3.08 MHz sampled with 63 channels for each 
IF (48.83 kHz, or 0.7 km s~^ spectral resolution). Two IFs were used to observe right and 
left circular polarizations (RCP and LCP) centering the bandwidth at vlsr = ~60 km s^^, 
while the other two IFs (RCP and LCP) were centered at vlsr = +89.6 km s~^. Based on 
the results of the first epoch, we set the phase center of these observations at R.A.(J2000) = 
16'*58™27.3*, Dec.(J2000) = — 30°55'08", to include all observed maser components within 
1" from this center. 

For both epochs, source J1331+305 was used as the flux calibrator (assumed flux density 
~ 2.54 Jy), J1229+090 was the bandpass calibrator (bootstrapped flux density 20.5 + 1.1, 
26.4 ± 0.4, and 25.2 ± 0.4 Jy on 2006 March 22, 2006 March 26, and 2007 December 24, 
respectively) , while J1626-298 was the phase calibrator (bootstrapped flux density ~ 1.00 + 
0.09, 1.72+0.03, and 2.30+0.04 Jy on 2006 March 22, 2006 March 26, and 2007 December 24, 
respectively). Maps were obtained with robust weighting of visibilities (Robust parameter 
= 0) and deconvolved using the CLEAN algorithm. The resulting synthesized beams were 
~ 0.20" X 0.07" and ~ 0.7" x 0.25" for the first (A configuration) and second (B configuration) 
epochs, respectively, with their major axes oriented nearly north-south. 

Data were self-calibrated using a strong maser feature, and applying the obtained phase 
and amplitude correction to the rest of the channels for each data set. In the first epoch, the 
maser components detected in each day (2006 March 22 and 26, i.e., blue and redshifted emis- 
sion, respectively) could be aligned to a common reference position, using spectral features 
that were present in the overlapping region between each individual observation. However, 
the data taken in different days could not be aligned since there was no common spectral 
feature in the velocity ranges covered. Therefore, the 2a relative positional uncertainty be- 
tween the maser features observed on the same day in the first epoch is < 0.001" (estimated 
as 2a ~ where is the size of the primary beam and SNR is the signal to noise ratio), 
while the positional uncertainty among the emission observed in different days is ~ 0.1". For 
the second epoch, since all data were taken simultaneously, the self-calibration corrections 
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obtained for the strongest maser feature was applied to all channels. The relative positional 
accuracy among all features in that second epoch was < 0.005" 

3. Results 

The spectrum and the spatial distribution of the water maser components obtained on 
the first epoch are shown in Fig. [H By "maser components" we mean individual intensity 
peaks in the spectrum. The positional information has been obtained by fitting elliptical 
gaussians to the maser emission, only in the channels in which a spectral peak is present. 
The total velocity span of the maser emission is ~170 km s"^. Two groups of components are 
evident in the spectrum: the redshifted group, from 50 to 110 km s~^, and the blueshifted 
one, from —40 to —80 km s~^. The red components are more numerous and more intense 
than the blue ones, with the most intense peak reaching ~ 5.1 Jy. The two groups appear 
spatially separated by ~ 0.11", forming a bipolar distribution, with the red components 
clustered at the NE and the blue components clustered at the SW. However, as explained 
in Sec. |2l while the spatial distribution of components within either group is very reliable 
(relative positional errors < 0.001"), the relative position between both groups has an error 
of ~ 0.1" (since they were observed in different days), which is on the order of the apparent 
separation. Therefore, the bipolar distribution of the maser emission observed in the first 
epoch may not have been real. 

Fig. [2] shows the spectrum and the spatial distribution of the water maser emission 
obtained in the second epoch. These observations did not cover the whole velocity range in 
which emission is present. They were set up to simultaneously observe red- and blueshifted 
components, in order to confirm the bipolarity. In these data, the NE-SW bipolar distri- 
bution is present, with a separation of ~ 0.08" between the centroids of the red and blue 
groups along position angle (PA) ~ 44°. The separation and PA are both more reliable for 
this second epoch (relative position error between features ^ 0.005"), although the internal 
distribution of the components in either group is more accurate in the first epoch. Therefore, 
water maser emission from IRAS16552 seems to trace a bipolar outflow with a projected size 
of ~ 0.08", and velocity along the line of sight of ~ 85 km s~^ (half of the total velocity 
span). Our results confirm IRAS16552 as a water fountain with a bipolar distribution of its 
water maser emission. It also confirms the trend of bipolarity in all water fountains studied 
up to now at high angular resolution with radio interferometry, strongly suggesting that 
water masers in all water fountains are produced in coUimated, bipolar jets. 

The centroid position of the individual water maser components has absolute coordinates 
R.A.(J2000) = 16^58'"27.30^ Dec(J2000) = -30°55'8.0" (absolute position error ^ 0.1", 
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compatible in both epochs). Assuming a bipolar configuration for the jet, we expect the 
powering source to be located near that position. Howe ver, this centroid is^ 16" away from 



the post-AGB star identified at optical wavelengths by lSuarez et al.l (120061 ) . indicating that 



the water maser emission is associated with a different evolved star within the field. Note 



that there is an error in the possition given in Table 3 of ISuarez et al.l (120061 ). since this 
position does not correspond to the KOI star marked in the identification chart and whose 
spectrum is shown in that paper, but to a nearb y source ('2MASS J 16582776-3055062). 



The true coordinates for the KOI star identified in ISuarez et al.l (120061 ) are R.A.(J2000) = 
16'^58™28.53^ Dec(J2000) = -30°55'09.3". 

We have searched the DSS plates and the 2MASS catalog for an optical and/or near 
IR counterpart at the position of the outfiow source, but although several optical and near- 
IR sources lie within the IRAS ellipse error, none of them is coincident with the position 
of the maser emission. The closest 2MASS source, J16582725-3055108, does not coincide 
within the errors (0.5" for 2MASS) with the maser position, since it is located at 2.3" from 
the centroid position. This implies that the outfiow source must be heavily obscured by 
its circumstellar envelope. From these results, it is plausible that the IRAS point source 
actually traces this obscured source that is powering the masers, but it is not related to the 
optical source. In this case, the IRAS colors will classify IRAS16552 as a post-AGB star, a 
classification also supported by its low IRAS variability. Its likely post-AGB nature, together 
with its association with high-velocity bipolar water-maser emission further indicates that 
IRAS 16552 belongs to the water fountain class. 

As for the spatial distribution of water masers within each group, the redshifted (NE) 
cluster of masers is oriented almost along the north-south direction and extends ^ 0.04"- 
0.05" in both epochs. The SW group is also oriented almost north-south in the first epoch 
and extends ~ 0.03". The distribution in the second epoch looks somewhat different, with 
an orientation along PA ~ 45°, but given the relative errors (see Sec. [2]), the distribution in 
this second epoch is less reliable. In any case, it seems clear that, if we assume that the line 
joining the red- and blueshifted groups traces the jet direction, the maser emission within 
each group is not elongated along this jet. 

As shown in Figs. [T] and [H the maser components within the redshifted cluster tend to be 
grouped according to their velocities, although the apparent grouping seems to be different in 
the two epochs. In the first epoch, the least redshifted group, with v < 60 km s~^ is located at 
the north, almost all the components of the intermediate- velocity group (60 < v < 80 km s~^) 
are located at the south, and the most redshifted components {v > 80 km s~^) are situated 
in between the other two groups. In the second epoch, the group centered ~ 75 km s~^ is 
north of the one centered ~ 90 km s~^, which is the opposite distribution to the one in the 
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first epoch. We note again that in the second epoch we did not observe the same velocity 
range as in the first one, and the positional error among the components of each group is 
larger in the second one (~ 0.005" as compared to ~ 0.001"). In the SW group, no clear 
tendency in the distribution of radial velocities is seen in either epoch. 



Discussion 



4.1. A comparison of IRAS16552 with other water fountains 

The absence of ne ar IR counterpart in the Two Micron All Sky Survey (2MASS) catalog 
( jSkrutskie et al.ll2006l ) does not seem to be an uncommon cha racteristic i n the water foun- 



tains. Of the 11 water fountains reported so far (tho se listed in 



by the Midcourse Space Experiment (MSX) catalog (lEgan et al 



mai 



20071 ) all those covered 



20031 ) (only IRAS16552 was 



not included in the MSX survey) have a mid-infrared counterpart. However, we have found 
2MASS counterparts (coincident with MSX sources within the errors) for only 5 of them. We 
did not find any significant relationship between the presence of a near IR counterpart and 
the characteristics of the water masers, such as number of maser components, dynamical age, 
or velocity span. The still small number of water fountains known is an obvious hindrance 
for this type of comparison. However, the fact that about half of the water fountains are not 
detected in the near IR gives us hints about the characteristics of their envelopes compared 
to those of the rest of AGB and post-AGB stars. Both AGB and post-AGB stars are usually 
detected in the near IR, alth ough a large fractio n of AGB stars, and some post-AGB stars are 
not detected in the optical (ISuarez et al.l 120061 ). This suggest that the central stars of water 
fountains are surrounded by thicker envelopes than in the rest of AGB and post-AGB stars. 
This could be explained if the evolution into the post-AGB stage is faster in water fountain 
sources, and therefore, these stars are relatively massive (^4 — 5 Mq). Nevertheless, a study 
of the near IR properties of water fountains, with a higher sensitivity than the data in the 
2MASS catalog may yield the detection of weak counterparts, provided that an accurate 
position of the water maser emission is determined with radio interferometric observations. 

In order to place IRAS16552 in the evolutionary framework of water fountains, we 
have plotted all such sources known on an IRAS two-color diagram (Fig. [3D, where we also 
show their dynamical ages (for those known). It is interesting that water fountains tend 
to be distributed i n this IRAS color-color diagram on the same re gion as post-AGB stars 



([Suarez et al.ll2006l) . while in this region, few AGB stars are located (jJimenez-Esteban et al. 



20061 ). This is certainly reasonable, given that water fountains seem to be post-AGB or 
evolved AGB stars. However, it is potentially important, since it suggests that IRAS colors 
could be useful as diagnostics to select potential candidate sources in future searches for 
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water fountains among the most extreme AGB stars. 



To try to establish the evolutionary stage for each of t he water fountai ns we have used 
the [8]-[12] vs [15]-[21] MSX two-color diagram (see FigS]). ISevensterl ^QM) has studied the 
position of OH-maser-emitting AGB and post- AGB sources in this diagram. The diagram 
has been divided in four quadrants (QI, QII, QUI, QIV), with classical AGB stars located 
in QUI (low left), early post- AGB stars in the QIV (low right) and more evolved objects in 
the QI, while QII contains star forming regions. The two water fountains with the shortest 
dynamical age are located in region III, and the one with the longest in region I, which 
agrees with Sevenster's evolutionary distribution. Estimates of the dynamical ages of the 
rest of water fountains will be helpful to confirm this trend. Moreover, the objects with 
visible lobes (IRAS 16342-3814 and I RAS 19134 - f-2131 ) are in the QIV and QII region, 
always with [8] — [12]>0.9. According to ISevensterl (120021 ). the transition from [8] — [12]<0.9 
to [8] — [12]>0.9 implies the transition from AGB to post-AGB stars. 

The water maser components in water fountains seem to show always a spatial bipolar 
layout, although their distribution is rather different from one object to another. Also the 
number of components is very variable from one source to another. There are sources as 
OH12.8— 0.9 which show only 2 or 3 water maser corn ponents with very similar velocities 
at the tips of each lobe jBoboltz fc Marvel! boosl . l2007h . and sources as IRAS 18286-0959 
that show multiple components with significant velocity dispersion within each lobe, and also 
components located close to the center of the maser distributio n, with intermediate velocities 
between those of the most extreme components (jimail 120071 ). In the case of IRAS16552, 
it shows at least 15 components with a spread in velocities within each lobe of ~ 25 — 
50 km s~^, and an absence of equatorial masers. These spectr al and spatial char acteristics 
make it simil ar to other water fount ains as IRAS 16342—3814 (IMorris et al.ll2003l ). or IRAS 
19134+2131 Jimai et al.ll2"ooi . l2007h . 



We note that single-dish obser vations of the OH maser line at 1612 MHz did not show 
any maser emission in IRAS16552 (ITe Lintel Hekkert et al.lll99ll : iHu et al.lll994l ). with la 
rms of 0.2, and 0.1 Jy, respectively. The presence or not of OH maser emission, as well 
as its kinematical and spatial characteristics is of potential interest to further characterize 
water fountains. The presence of OH maser emission is a common characteristic in almost all 
water fountains. Apart from IRAS16552, the o nly other known water fountain iri which OH 
seems not to be present is IRAS 19134-1-2131 Jhewis et al.lll987l : IHu et al.lEggl . The OH 
masers detected are in general less extended and have a smaller range of velocities than the 



H2O masers, as is th e case of IRAS 18043—21 26 (ISevenster fc Chapmanll200ll : iDeacon et al. 



2007h . OH12.8-0.9 JBoboltz fc Marvell l2005h . W43A Jimai et al.ll2002h . This suggest that 
the OH masers tend to trace the former AGB shell, rather than the later bipolar mass-loss 
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traced by H ^O. However, the cases of IRAS 15445-5449 and IRAS 15544-5332 f lDeacon et al. 



2004J . 120071 ) do not follow this trend, since OH maser emission in these sources have been 



found outside the velocity range of the H2O masers. IRAS 16342—3814 can be considered 
as another exception since although t he OH masers have a smaller velocity range than the 
H2O masers (70 krn s~^ y s 160 km s~^; ISahai et al.lll999l : IZijlstra et al.ll200ll : lLikkel fc Morris 



19881 : iMorris et al.ll2003l ). they are located in the bipolar jets and have a velocity that does 



not correspond to the AGB shell. 



4.2. The orientation of the water maser distribution 



As we mentioned in Sec. |3l the water maser emission in IRAS16552 is not aligned 
with the apparent iet direction. Such mi salignment is also common in other sources, as 
IRAS 19134+2131 fllmai et all booi l2007h and W43A Jlmai et al. I2OO2I). Those cases ca n 
be explained as the effect of ballistic corkscrew jets (jimail 120071 : IVlemmings et al.ll2006l ). 
where the masers are emitted in a ballistic way along the direction of a jet that precesses 
with time. This precession could be due to the presence of a binary stellar companion or 
a massive planet. These companions wo uld also account for the magnetic fields found in 
sources as W43A (jVlemmings et al.ll2006l ). Nevertheless, no direct evidence for a companion 
has been found so far in any water fountain. 

IRAS16552 is an extreme case in which the water maser components are oriented almost 
perpendicular to the radial direction from the central source (assumed to be at the center of 
the distribution of the comp onents, see Sec. [3D . An almost perpen dicular distributiori is als o 
seen in IRAS 16342-3814 jMorris et al.ll2003h and OH 12.8-0.9 JSoboltz fc Marvel! l2005h . 
which has been attributed to the masers tracing a shock front. For IRAS16552, such an 
interpretation would imply a rather large opening angle of the jet (~ 70°). A possible alter- 
native explanation is that the water masers in IRAS16552 trace a rapidly precessing/rotating 
bipolar jet. If the precession/rotation period is very short, as compared with the expansion 
timescale determined by the jet velocity, the bipolar jet (or the shocked regions in the neu- 
tral envelope) will describe point-symmetric structures oriented almost perpendicular to the 
outflow direction. At least qualitatively, this scenario may account for the observations of 
IRAS16552, specially for the first epoch, in which the distribution of masers within each 
group has been determined with a higher accuracy. 

Remarkably, many PNe show coUimated structures oriented almost perpendicular to 
the radial direction from the central star. This is the case of bipolar and elliptical PNe 
that exhibit a poiiit-syra metric inten sity distribution in the shell as, for instance, Hb5 
( ICorradi &: Schwarzlll993l ) and Cn3-1 (IMiranda et al.l 119971 ). More interesting is the case of 
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components DD' in NGC6543 (IBalick fc Hajianll2004l ). which are dense (10^ - 10^ cm~^) fila- 
ment s that have been interpreted as due to a rapidly precessing bipolar jet (IMiranda fc Soli 



19921 ). Given that water fountains are probably related to relatively massive stars (Sec. 14.1 



their evolution through the post-AGB phase should be very fast. The fingerprints of rapidly 
precessing, collimated post-AGB outfiows, such as the one proposed for IRAS16552, may 
still be present in the envelope when the central star ionizes the nebula and could show up 
as point-symmetric ionized filaments in PNe, similar to the components DD' in NGC6543. 



5. Conclusions 

In this paper we have presented VLA observations of H2O masers at 22 GHz towards 
IRAS 16552—3050. Our main conclusions are as follow: 

• We have found that the water maser emission in IRAS 16552—3050 span a velocity 
range of ~170 km s~^ and shows a bipolar distribution of ~ 0.08" in size. The properties 
of the water maser emission and the likely post-AGB nature of the central source allow 
us to confirm IRAS 16552—3050 as a member of the class of "water fountain" evolved 
objects. 

• Our result confirms the trend of bipolarity in all water fountains known up to now, 
strongly suggesting that water masers in these sources are produced in bipolar, colli- 
mated jets. 

• The water maser emission in I RAS 16552—3050 doe s not seem to be ass ociated with 
the optical source identified by iHu et al.l (119931 ) and ISuarez et al. J2006h . but with a 
different evolved object, for which there is no optical nor near-IR counterpart known. 

• About half of the water fountains known do not have a near-IR counterpart in the 
2MASS catalog. This suggest that water fountains are surrounded by thicker envelopes 
than the rest of AGB and post-AGB stars (which usually have near-IR counterparts), 
probably implying that they are relatively massive objects (^4 — 5 Mq). 

• The distribution of the water maser emission in the red- and blueshifted group is almost 
perpendicular to the proposed jet direction. This might be explained if the maser 
emission traces a rapidly precessing/rotating bipolar jet. A later evolution could give 
rise to point-symmetric filaments observed in some PNe, which appear oriented almost 
perpendicular to the radial direction from the central star. 
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c) 




Fig. 1. — (a) Spectrum of the water masers towards IRAS 16552—3050, observed on 2006 
March. The different maser components are identified by colored stars. The color code is 
used to identify the same components in the rest of the panels, (b) Close up of the redshifted 
part of the spectrum, (c) Spatial distribution of the peaks of the maser components. Note 
that the coordinates in the figure should be considered only as nominal. The error in absolute 
position is ~ 0.1" (see Sec. [2]). (d) Close up of the spatial distribution, showing only the 
redshifted components. The symbol size increases with flux density. 
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Fig. 2. — a) Spectrum of the water masers towards IRAS 16552—3050, observed on 2007 
December. The different maser components are identified by colored stars. The color code 
is used to identify the same components in panel b. (b) Spatial distribution of the peaks of 
the maser components. Note that the coordinates in the figure should be considered only 
as nominal. The error in absolute position is ~ 0.1" (see Sec. [2]). The symbol size increases 
with flux density. 
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Fig. 3, 

(see 



Positi on of the 11 water fountains known up to now in the IRAS two color diagram 



Imail 120071 . and references therein). The dynamical age is shown between parentheses 



for those objects for which it has been calculated. The curve in the plot is the one defined 
by iBedijnl ( 119871 ) where the AGB stars tend to be located. 
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Fig. 4. — MSX two-color diagram of the 10 water fountains known up to now tha t have been 
observed with this sateUite. The four quadrants are those defined in ISevensterl (120021 ) and 
described in the text. 



